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Abstract: We revisit ghost dark matter, the possibility that ghost condensa- 
tion may serve as an alternative to dark matter. In particular, we investigate the 
Friedmann- Robertson- Walker (FRW) background evolution and the large-scale struc- 
ture (LSS) in the AGDM universe, i.e. a late-time universe dominated by a cosmolog- 
ical constant and ghost dark matter. The FRW background of the AGDM universe 
is indistinguishable from that of the standard ACDM universe if M > 1 eV, where 
M is the scale of spontaneous Lorentz breaking. From the LSS we find a stronger 
bound: M > 10 eV. For smaller M, ghost dark matter would have non-negligible 
sound speed after the matter-radiation equality, and thus the matter power spectrum 
would significantly differ from observation. These bounds are compatible with the 
phenomenological upper bound M < 100 GeV known in the literature. 
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1. Introduction 

Current data of the cosmological observations (e.g., Cosmic Microwave Background 
(CMB), Large Scale Structure (LSS) and SuperNovae (SNe)) show that our Universe 
today mostly consists of dark matter responsible for the structure formation and dark 
energy causing late time accelerated expansion of the Universe |I] 0, |3], [|, [| . From 
observational point of view, the paradigm of dark energy and dark matter is very 
successful to fit the data. 

However, from theoretical viewpoint, we do not know what they really are, de- 
spite the fact that there are many theoretical models (for review, see e.g. [|] for 
dark energy and |7]] for dark matter). This situation has been a strong motivation 
for modification of gravity as an alternative to dark energy and dark matter: just 
changing behavior of gravity at long distance/time scales might be able to explain 
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the observational data without introducing dark energy and dark matter. Indeed, 
many theories of modification of gravity have been proposed, such as massive grav- 
ity p, DGP model || [T0| , ghost condensate [|11|, [12] and so on. It is important to 
investigate cosmological implications of those modified gravity theories toward the 
future observations. 

In this paper, we focus on the ghost condensate scenario and investigate the pos- 
sibility that it might serve as an alternative to dark matter. This possibility, dubbed 
ghost dark matter, was already pointed out in [|ITJ but has not been investigated 
in detail. Based on the Friedmann-Robertson- Walker (FRW) background evolution 
and the large-scale structure of the universe, in the present paper we shall find a 
lower bound on the scale of spontaneous Lorentz breaking, M > 10 eV, under the 
assumption that ghost dark matter is responsible for all dark matter in the universe. 
Most importantly, this bound is compatible with the phenomenological upper bound 



M < 100 GeV found in 12 



The rest of this paper is organized as follows. In the next section, we briefly re- 
view the ghost condensate scenario, including the basic idea, the low-energy effective 
theory and the phenomenological upper bound on the scale of spontaneous Lorentz 
breaking. In Sec. [|, we introduce a simplified description of ghost dark matter and 
investigate the FRW background evolution. We also clarify the regime of validity of 
the simplified description. In Sec. f|, we consider effects of ghost dark matter on the 
large-scale structure of the universe. We discuss how density perturbations evolve 
in the universe dominated by the cosmological constant and ghost dark matter, i.e. 
the AGDM universe, and then give a constraint on the model from the shape of the 
matter power spectrum. The final section is devoted to conclusion of this paper and 
discussions. 



2. Review of ghost condensation 
2.1 Basic idea 

In particle physics it is the so-called Higgs mechanism that modifies force law in 
the infrared (IR) and that makes it possible to describe the weak interaction in a 
theoretically controllable way. A non-vanishing vacuum expectation value (vev) of a 
scalar field spontaneously breaks a part of the gauge symmetry and modifies the IR 
behavior of the corresponding gauge force from Gauss law to Yukawa law. 

Ghost condensation applies the idea of Higgs mechanism to general relativity to 
modify the IR behavior of gravity [11], 0. In order to spontaneously break a part 
of the symmetry of general relativity, i.e. spacetime diffeomorphism invariance, we 
consider a non-vanishing vev of derivative of a scalar field. In addition, we demand 
that the vev is timelike so that only the time reparametrization symmetry is spon- 
taneously broken and the 3-dimensional spatial diffeomorphism invariance remains 
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unbroken. Note that this symmetry breaking pattern is totally consistent with obser- 
vational fact that our universe averaged over large scales is isotropic. However, unlike 
usual situations in cosmology with a scalar field, we require that the vev of derivative 
should remain non-zero and finite as the universe expands towards maximally sym- 
metric spacetimes, i.e. Minkowski or de Sitter spacetimes 1 . This is because we would 
like the IR modification of gravity to persist in Minkowski or de Sitter backgrounds. 

Therefore the equation of motion of the scalar field and the Einstein equation 
must allow a solution with X = —d^tfid^cj) constant and positive 2 in Minkowski or 
de Sitter spacetime. This requirement forbids inclusion of any non-trivial potentials 
for the scalar field since a potential for an eternally running scalar field would lead 
to time- dependence in the stress-energy tensor. In other words, we must invoke the 
shift symmetry for the scalar field action, i.e. invariance of the action under constant 
shift of the scalar field: — > <f) + c, where c is an arbitrary constant. With the shift 
symmetry, the action for the scalar field minimally coupled to gravity should be of 
the form 



where Y = VV^V^V^, Z = V^V" (f>d^(f)d v (f), and so on. In addition to the shift 
symmetry, we assume the Z2 symmetry for the scalar field action, i.e. invariance of 
the action under reflection <ft — > —(f). The Z 2 symmetry is to ensure that the effective 
theory for excitations of ghost condensate is invariant under simultaneous reflection 
of the time t and the scalar field perturbation 5(f): t — > —t, 5(f) — > —5(f). 
The simple Lagrangian 



i.e. either Px — > or X — > 0. Note that X is either positive or zero for a homogeneous 
(f). Thus, if P(X) has an extremum at X = M 4 > 0, i.e. P X {M A ) = 0, then 
X = M 4 is a dynamical attractor of the system and ghost condensate may be realized 
automatically. 

However, we shall see below that a more general action depending on e.g. □</> 
(see ( |2.1| )) is needed to describe excitations around the ghost condensate properly. 

1 We do not consider anti-de Sitter spacetimes as backgrounds since we are interested in cosmol- 
ogy. If we nonetheless considered ghost condensation in an anti-de Sitter background then the vev 
of derivative would be spacelike and excitations around the condensate would not have a healthy 
kinetic term. 

2 In this paper we adopt the mostly positive sign for the spacetime metric. 
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2.2 Decoupling limit in Minkowski background 

Let us now consider excitations around the extremum X = M A > of P(X). For 
simplicity we consider Minkowski background. We shall expand the action (|2.2j ) 
around (j) — ct and take the limit c 2 — > M 4 . Note that in Minkowski background the 
equation of motion for is satisfied for any c and that this treatment is consistent 
as far as the backreaction to the geometry is negligible. 
For 

= ct + ir(t,x), (2.5) 
by expanding P(X) with respect to 7r we obtain the quadratic Lagrangian for 7r as 

Lf = [2c 2 P xx (c 2 ) + P x (c 2 )] (d t n) 2 - P x (c 2 )(Vrc) 2 , (2.6) 

where Pxx = d 2 P/dX 2 . Note that we did not take into account the backreaction of 
the scalar field to the background geometry nor include metric perturbations. These 
treatments are justified in the limit E/Mp\ — > 0, where E represents energy scales 
of interest and Mpi is the Planck scale. Since we shall later take the limit c 2 — > M 4 
and we shall see that M sets the cutoff scale of the effective field theory, E can 
be replaced by M in the regime of validity of the effective field theory. Thus, the 
decoupling limit is characterized by M/Mp\ — > 0. In this limit and for energies and 
momenta sufficiently lower than M, the action (2^6) is valid and the small fluctuation 
7r is stable if 

2c 2 P xx (c 2 ) + P x (c 2 ) > 0, P x (c 2 ) > 0. (2.7) 
By taking the limit c 2 — > M 4 , we obtain 

L^ = 2M 4 P xx (M 4 )(d t n) 2 . (2.8) 

Thus, the excitation it around the attractor X = M 4 has a healthy time kinetic term 
if P XX (M A ) > 0, i.e. if X = M 4 is a local minimum of the function P{X). However, 
the spatial gradient term vanishes. 

This means that the action ( |2.2|) is too simple to describe excitations around the 
ghost condensate background. Indeed, while we have assumed the shift symmetry 
to prevent a non-trivial potential from being generated, there is no way to prevent 
□0, Y, Z, etc. from appearing in the action. We therefore have to go back to the 
general action (2~T) and seek leading gradient terms. 

It turns out that the leading gradient term is of the form 

A^ = -^(V 2 vr) 2 (2.9) 

where a is a constant of order unity. (The reason why this is indeed the leading term 
will be made clear in the next subsection.) Combining this with the leading time 
kinetic term ( |2.8|) , we obtain the quadratic action for 7r 

dtd 3 xL n , = M A ^-{d t ix) 2 - ^ 2 {V\) 2 V (2.10) 
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where we have normalized it and a as 2^/ ' Pxx(M A )ir — > ir and a — > AM 4 Pxx(M 4 )a. 
Thus, the dispersion relation for it in the decoupling limit M/Mpi — > is 



«■ = (2.11) 



2.3 Scaling dimension and suppression of extra terms 



The general action ( |2.1|) should in principle include any terms consistent with the 
shift symmetry and the Z2 symmetry. In terms of n, they are invariance under the 
constant shift of 71 (tt — > it + c, where c is an arbitrary constant) and invariance under 
the simultaneous reflection of the time t and it (t — > —t, n — > — it), respectively. Thus, 
there should be infinite number of terms added to the quadratic action ( |2.1(J| ) for n. 
Nonetheless, one can show that the quadratic action fl2.10| ) is a good description of 
low energy behavior of the system. 

In order to show that those extra terms are irrelevant at low energies, let us first 
identify the scaling dimensions as the energy E is scaled by E — > sE (or the time 
interval dt is scaled by dt — > s _1 c/t), where s is some constant. By requiring that the 
quadratic action ( [2.10| ) be invariant under scaling, we can fix the scaling dimensions 
as 

E ->• sE, 
dt — > s _1 <it, 
dx — > s~ 1 / 2 dx, 

7T ->• S 1/4 7T. (2.12) 

Note that the scaling dimensions of dt and dx are consistent with the dispersion 
relation ( |2.1ip . With this scaling, one can check that the leading interaction 



M 4 J dtd 3 x(V7i) 2 d t 7i (2.13) 

scales as s 1//4 . Thus, this term is irrelevant and becomes less and less important 
at energies and momenta sufficiently lower than M. All other terms are even more 
irrelevant. 

There is one relevant operator, namely (V7r) 2 . However, as we have already seen, 
the coefficient of this operator is proportional to Px and goes to zero as the universe 
expands. 

Therefore, if energies, momenta and the field amplitude are sufficiently lower than 
M then low energy/momentum/amplitude expansion around the quadratic action 
( p. 10 ) is under control. This in particular implies that apparent extra modes due to 



higher time derivative terms have frequencies of order M or higher. Hence there is 
no ghost in the regime of validity of the effective field theory if we set the cutoff scale 
slightly below M. 
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2.4 Jeans-like instability and IR modification of linear gravity 

So far, we have considered the decoupling limit M/Mp\ — > of the theory. For small 
but finite M/Mp\, the dispersion relation fl2.11|) gets corrected due to mixing with 
gravity and becomes [[ll| 

This dispersion relation exhibits Jeans-like instability for modes with length scales 
longer than L c and the corresponding time scale is T c , where 

L.„% T c ~-Mk-. (2.15) 

Note that this is an IR instability and has nothing to do with ghost. The Jeans- 
like instability is the origin of the IR modification of gravity in ghost condensate 
background [pTlJ] . The time scale and the length scale of the modification are T c and 
L c , respectively, and are much longer than the naive scale 1/M. 

In Sec. |4| we shall investigate Jeans instability of ghost dark matter, a component 
which arises from excitation around ghost condensate and which behaves like dark 
matter. One should note that the Jeans-like instability in the exact ghost conden- 
sate background considered in this subsection is both conceptually and qualitatively 
different from the Jeans instability of ghost dark matter in Sec. |j. 

In the decoupling limit M/Mp\ — > 0, the timescale T c diverges (in the unit of 
1/M) and thus the Jeans-like instability disappears. 

If we required that the Jeans timescale T c be longer than the age of the universe 
then we would end up with the (would-be) upper bound M < 10 MeV JTT]. However, 



we shall see below that nonlinear dynamics becomes important much earlier than T c , 
that this (would-be) bound is not necessary and that the current upper bound on M 
is as weak as M < 100 GeV |2j. 

2.5 Nonlinear dynamics and upper bound on M 

A sightly nonlinear extension of the quadratic action ( |2.10|) coupled to linearized 
gravity is 



dtd 3 xL n , L n = M 



2 



0,K-(YKf-i> -^(VV!>. (2.16) 



where $ = —Sgoo/2 is the Newtonian potential |12| in the longitudinal gauge. This 



is obtained by not dropping the leading nonlinear term and the metric perturbation 
in going from ( |2.2| ) to (|2.8|) before adding (|2.9|). Note that (d t 7i) 2 has been dropped 
from the expression in the square bracket since it has a scaling dimension higher 
than (V7r) 2 , following the discussions in subsection 
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From the quadratic part of the action ( |2.16|) , i.e. (|2.10|) , it is easy to see that 
the timescale of linear dynamics Tun is determined by 



7T 2 7T 2 



(2.17) 



where we have assumed that the length scale of interest L is shorter than the Jeans 
scale Lj (so that the second term in the right hand side of ( |2.14| ) is negligible) and 
we have set a = 0(1). Thus we obtain 

T Lin ~ ML 2 , (2.18) 

which is consistent with the dispersion relation fl2.11|) . On the other hand, from the 
nonlinear action ( 2.16Q , the timescale of nonlinear dynamics Tnl is determined by 



7T 7T 2 



and we obtain 



$, (2.19) 



MIL 3 

T ^~^wri^t- (2 ' 20) 

where M src is the mass of the gravitational source. Note that this timescale is noth- 
ing but the Kepler time. Therefore, nonlinear dynamics dominates before linear 
dynamics if T NL < T Lin , i.e. 

m ^ I r 

> W^^. 2.21 



M P1 ~ V M BIC L 

This condition is satisfied in virtually all interesting situations. For example, for the 
earth's surface gravity, this condition is as weak as 

M > 1(T 9 eV. (2.22) 

From the dispersion relation ( |2.14 ) it is easy to see that modes with L > Lj 
grow due to Jeans-like instability in timescale 

r ~ ^L. (2.23, 

The nonlinear term (Vvt) 2 in the squared bracket in the action (|2.16|) becomes im- 
portant when it is comparable to <9j7r or larger, i.e. when 

|tt| > tt c = — . (2.24) 
r 

This can be rewritten as a condition on the energy density p n ~ M 4 d t 7r: 

__M^ C M 6 
P ^ Pc = — ~Ml{ (2 ' 25) 
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Hereafter, we assume that nonlinear dynamics cutoff the Jeans-like instability at 
|Ptt| ~ Pc- Positive and negative regions can grow to \p n \ ~ p c within the age of the 
universe Hq 1 if L > L c and r < Hq 1 , i.e. if 

L c < L< L max = — ~ 10R G x — — — , (2.26) 
c ~ ~ max M Pl H ° lOOGeV V ; 

where R Q is the solar radius. Therefore, if M is lower than 100 GeV (we shall indeed 
see below that M must be lower than 100 GeV) then the typical size of positive and 
negative regions is smaller than 10R & . This means that p n averaged over large scales 
such as galactic scales is almost zero and does not gravitate significantly. 

Based on these properties of nonlinear dynamics, several phenomenological upper 
bounds on M were derived in |T2| . The strongest among them is the 'twinkling- from- 
lensing' bound, which we shall briefly describe here. 

Suppose that the universe is filled with regions with p^ ~ ±p c of the size L c < 
L < L max moving relative to the cosmic microwave background (CMB). Each region 
contributes to weak gravitational lensing with the deflection angle 



r„ M% 2 M 6 L 2 

Atfcach ~ — ~ -a a-, 2.27 

Mpj Mpj 



where b (< L) is the impact parameter and r g ~ p c b 3 /Mp X is the gravitational radius 
of the mass contained within the impact parameter. In the final expression, we have 
maximized A# cach with respect to b. Since a light-ray from the distance d experiences 
N (~ d/L) lens events, the total deflection angle is 

AO AO W/2L3/2 W/2L - (22R] 
A0 tot ~ A0 cach VN ^ ^ , (2.28) 



where we have maximized A0 tot with respect to L in the final expression. We can 

to 



apply this result to the CMB by setting d ~ H . Requiring that A0 tot be smaller 



than the angular resolution of CMB experiments ~ 10 3 , we obtain the upper bound 

M < 100 GeV. (2.29) 

The twinkling timescale for the CMB is 

L max M 300km/s 

twinkle innn „ — x 0.1 day, 2.30 

v lOOGeV v 

where v is the typical velocity of positive and negative regions relative to the CMB 
rest frame. Thus, if M is close to 100 GeV then the twinkling effect may be detected 
in future CMB experiments. 
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3. Ghost dark matter and background evolution 



3.1 Simple description of ghost dark matter 



Now let us reconsider the simple Lagrangian (|2.2| ) which depends on X only. As 
we have already seen in subsection |2.2| , we need to add the extra term ( |2.9|) to this 
Lagrangian in order to describe perturbations around the exact ghost condensate 
background, i.e. a local minimum of the function P{X). In this subsection we shall 
instead consider a background with non- vanishing Px and see that, under a certain 
condition, the Lagrangian P(X) without the additional term ( |2.9| ) can properly de- 
scribe the background and perturbations around it. In this situation, as we shall see 
below, deviation from the exact ghost condensate behaves like dark matter. 
The stress-energy tensor corresponding to the Lagrangian P{X) is 

T„ v = 2P x d lM (j)d v (j) + P 9fMU = (p + P) Uy u v + P 9lu/ , (3.1) 

where 

p = 2XP x -P, Ua4 = _M. (3.2) 

v X 



The sound speed squared for perturbation is Jl3| 



= dP/dX = P x 
s dp/dX 2XP XX + Px' ' 

This agrees with minus the ratio of the coefficient of the gradient term to the coeffi- 
cient of the time kinetic term in the quadratic Lagrangian ( |2.6|) . 

As we have already seen in subsection [2.2| , in order to describe perturbations 
around the exact ghost condensate background, the simple Lagrangian P(X) is not 
sufficient but we need to add the term ( |2.9[ ) to it. The reason is that the coefficient 
of the would-be leading gradient term ( W) 2 vanishes in the exact ghost condensate 
background. However, if the background is not exact ghost condensate but has 
Px 7^ then the coefficient of (W) 2 does not vanish as shown in Q3.3|) . Therefore, 
if c 2 is large enough then we do not have to add the extra term ( |2.9|) to the simple 
Lagrangian P(X). To be more precise, the extra term (|2.9|) is not needed if 



< » mw> (3 ' 4) 

where L is the length scale of interest and we have supposed that a = 0(1). 

We have also seen that a local minimum of P{X) is a dynamical attractor in 
the expanding universe. Thus, it is rather natural to Taylor expand P(X) around a 
local minimum X = M 4 as 

P ~ P(M 4 ) + ^P XX (M 4 )(X - M 4 ) 2 . (3.5) 
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This expansion of P{X) is valid if and only if \X — M 4 | <C M 4 . Note that we any- 
way have to restrict our consideration to this regime; otherwise, higher dimensional 
operators, which in general depend on not only X but also □</>, Y, Z, etc., would 
be unsuppressed and the system would exit the regime of validity of the low energy 
effective theory. (See discussions in subsection |2l|) Correspondingly, we have the 
following expansions. 



2 X-M 4 
s ~ 2M 4 

These expressions are rewritten as 



p ~ -P(M 4 ) + 2M 4 P XX {M 4 ){X - M 4 ), (3.6) 

(3.7) 



where 



P - Pgdc + Pgdm, P ^ Pgdc + Pgdm, C\ ~ (3.8) 
Pgdc = "Pgdc = COnst., Pgdm = -~. (3.9) 

Here, M 4 = 4M 8 P XX (M 4 ) ~ M 4 . In the regime \X — M 4 | <§C M 4 , it is intriguing 
to note that (p g d e , -Pgde) and (p g dm, -Pgdm) behave like dark energy and dark matter, 
respectively. In particular, we shall call the latter ghost dark matter [|12] . 

In terms of ghost dark matter component, the necessary condition ( |3.4| ) for the 
validity of the simple Lagrangian P(X) is written as 

M«L^.10-GeVx( I A-).(^) 1/2 , (3,0, 

where po = SM^Hq is the critical density. For example, if M 10 15 GeV and if we 
suppose that the ghost dark matter is responsible for all dark matter in the universe 
(p gdm ~ 0.3 x po) then (|3.10|) is satisfied for length scales longer than ~ 1 Mpc. 
Note, however, that validity of low energy effective theory requires that all associated 
energies, momenta and amplitudes be sufficiently lower than unity in the unit of M. 
Thus, ( |3.8| ) and ( |3.9| ) are valid only if 

Pgdm < M\ H < M. (3.11) 

On the other hand, we do not have to require that the radiation temperature be 
lower than M since interactions between ghost condensate and radiation are highly 
suppressed (typically by the Planck scale). 

3.2 Ghost dark matter production from ghost inflation 

In the previous subsection we have seen that the background evolution and the be- 
havior of perturbations of ghost dark matter can be described by (31!) and 
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under the conditions ( 3.10|) and ( p.llj ), where L is the length scale for perturba- 



tions of interest and M ~ M. If ( |3.10|) is not satisfied then the use of the simple 
Lagrangian of the form P(X) is not justified and we need to take into account ef- 
fects of the extra term ( |2.9|) . This is not a big problem but would make analysis 
slightly complicated Fortunately, for the purpose of the present paper, i.e. for 



understanding of the evolution of the FRW background and the large-scale structure 
of the universe, we are interested in L of order 1 Mpc or longer. In this case, the 
condition ( |3.10D always holds if the phenomenological upper bound M < 100 GeV 
(see subsection p75|) is satisfied and if we suppose that a non-trivial fraction of dark 
matter of the universe is ghost dark matter. On the other hand, if ( |3.11|) is not met 



then the system exits the regime of validity of the low energy effective theory and 
we need a ultraviolet (UV) completion 3 to describe the system. 

This naturally leads to the question "what happens in the early universe?" The 
condition (|3.11| ) does not hold in the very early epoch of the universe. In this early 
epoch the sector including ghost condensation should be governed by a theory more 
fundamental than what we have been describing since M is the energy scale above 
which a new physics kicks in. While it is important to seek a UV completion to 
describe this epoch properly, it is also plausible to consider cosmological scenarios 
in which all interesting observables are predicted within the regime of validity of the 
low energy effective theory. As a possible realization of such scenarios let us consider 
a generation mechanism of ghost dark matter at the end of ghost inflation |2T| . 



In ghost inflation the scalar field <fi responsible for ghost condensation plays 
the role of inflaton as well. For example we can consider a hybrid inflation-type 
implementation. In this case we suppose that the mass squared m 2 x of a water-fall 



field x depends on <p in such a way that m 2 ((p) is positive and constant for <p <C — 



'*3 



and negative and constant for (f> ^> 0*, where is a positive constant. This setup 
is technically natural since the shift symmetry is broken only in the vicinity of the 
transition region |0| < 0* and otherwise exact. We suppose that d t <fi > so that the 
sign of m? changes from positive to negative. 

A 

For (f) <C —(f)*, 4> en j°y s the shift symmetry and \ has a positive mass squared. 
Thus the system has a de Sitter attractor with d t <f) = M 2 and x — 0, where X = M 4 
is a local minimum of P(X). We suppose that the system settles in this state well 
before <fi crosses the transition region. Noting that the fluctuation S(p of has the 
scaling dimension 1/4 (see subsection |2.3| ) and the mass dimension 1, we can easily 
estimate the amplitude of quantum fluctuations of as 

6<P ~ M {^f) ' (3 ' 12) 



3 See (T^, [I?], |l8| for some attempts towards possible scenarios of UV completion, and [[l9| |2t 
for compatibility with the generalized second law of black hole thermodynamics. 
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where H m f is the Hubble expansion rate of the de Sitter attractor. As usual, quantum 
fluctuations of <fi is eventually converted to temperature anisotropics as 



ST H 



inft 



T d t c 

Thus, we obtain 



ST (H 



T V M 



Mnf 



(3.13) 



5/4 

(3.14) 



By requiring that this is responsible for the observed amplitude of temperature 
anisotropics ST /T ~ 10~ 5 , H in { is determined as 

H in{ ~ 1(T 4 x M. (3.15) 

Under this condition, one can also estimate non-Gaussian features of CMB anisotropies. 
The shape of the bispectrum is of the equilateral type and the nonlinear parameter 
J'nl is of order ~ 80 if we set all dimensionless parameters to unity |T!J. The es- 
sential reason for the relatively large non-Gaussianities is that the leading nonlinear 
operator has the scaling dimension 1/4 and thus is less suppressed than in usual 
slow-roll inflation. 



The condition ( |3.15| ) shows that ghost inflation is well within the regime of 



validity of the low energy effective theory: if inf <C M. The condition (|3.15|) also 



shows that there is a lower bound on M in terms of the reheating temperature T ie ^: 

T 2 

M > 10 4 x (3.16) 
Mpi ( J 

Ghost inflation can generate not only temperature anisotropies observed in the 
CMB but also ghost dark matter. In the hybrid inflation-type implementation, the 
shift symmetry is broken and a potential for <f> should be generated by quantum cor- 
rections in and only in the vicinity of the transition region \<f>\ < 0*. Correspondingly, 
the equation of motion for homogeneous is 

^d t [a^Pxdtcj)} + V g ' cn (0) = 0, (3.17) 

where V zen ((f>) is the potential for (p generated by quantum corrections. Note that 
it is not V gen ((j)) but a potential for \ that is responsible for most of the potential 
energy during the inflationary phase: 

lAVgJ < SM^H?^ (3.18) 

where AV gen = V^f ore — V^ cr , and the superscripts 'before' and 'after' represent 
values before and after the transition. Noting that d t <f) — M 2 arid ignoring the 
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Hubble friction term, it is easy to integrate ( 3.17 ) over the transition region. The 
result is 

P x teT - ^F- ( 3 - 19 ) 
This gives the amplitude of ghost dark matter right after the transition as 

~ A^gen- (3.20) 

Ghost dark matter does not dominate the universe just after the end of ghost inflation 



as easily seen from (|3.18| ). However, since it evolves more slowly than radiation, it 
can dominate the late time universe. 

In obtaining the estimate ( |3.19| ) we have ignored the Hubble friction. This is 
justified if the transition timescale, At ~ — <f)*/M 2 , is short compared with 

the Hubble timescale Hr* ~ 10 4 /M, i.e. if 

(j>* < 10 4 x M. (3.21) 

3.3 Phenomenological constraint from background evolution 

While it is certainly interesting and important to investigate concrete cosmological 
scenarios such as one presented in the previous subsection, it is also important to 
constrain the ghost dark matter in a model independent way. This subsection and 
the next section are devoted to this subject, assuming that ghost dark matter is 
responsible for all dark matter in the universe. 

We now consider background evolution from the radiation-matter equality to the 
present time to give a lower bound on M. In this epoch we know from observations 
that the background dark matter component behaves like pressure-less dust. Thus, 
if all dark matter in the universe is ghost dark matter then P g dm/Pgdm ~ Pgdm/M 4 
must be sufficiently lower than unity. (This condition incidentally agrees with the 
first inequality in ([3.11|) .) In this behaves as 

Pgdm - figdmPo • f — J • (3.22) 

We set f2 g dm — 0.3 by the assumption that ghost dark matter is responsible for all 
dark matter in the universe. Hence, by requiring that p g d m M A all the way up the 
matter-radiation equality a ~ 10~ 4 ao, we obtain the lower bound 

M > 1 eV. (3.23) 



4. Large-scale structure with ghost dark matter 

In this section, we calculate the evolution of density perturbations in the AGDM 
universe, where the sound speed for matter perturbation, c s , is given by Eq. (|3.7j ) 
unlike the standard ACDM universe. Then, we consider the bound on M in order 
not to conflict with current observations of large-scale structure. 
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4.1 Jeans wavenumber for ghost dark matter 



In subsection [3.3| , we obtained the lower bound on M as M > 1 eV by demanding 
that Pgdm/Pgdm "C 1 all the way up to the matter-radiation equality. The condition 



M > 1 eV also justifies the use of the simplified description given in subsection [Tl 
for the matter dominated era since -P g dm/Pgdm <C 1 is equivalent to p g d m *C M 4 . 

We consider the evolution of the matter perturbations in this situation. In this 
subsection, before detailed numerical analysis in the next subsection, let us foresee 
that a bound stronger than M > 1 eV can be obtained. 

For the sound speed given by Eq. (|3.7|), we can define the Jeans wavenumber as 



As is well-known, the matter perturbations do not evolve very much during radiation 
dominated era and we need not to consider the effects of the Jeans scales. However, 
after the matter-radiation equality we have to take into account the effects of the 
Jeans scales, namely, the perturbations with scales shorter than the Jeans scales 
can not grow. From Eq. ( |4.1| ), the Jeans wavenumber evolves as kj oc a during 
the matter dominated era and hence in the matter dominated era the comoving 
Jeans scale (~ kj 1 ) becomes the largest at the matter-radiation equality. The Jeans 
wavenumber at the matter-radiation equality is given by 

/f!, dm h 2 V 5/6 / M \ 4 ' 3 

The matter power spectrum is expected to be significantly suppressed for modes with 
wavenumbers k > kj teq . 

Almost all current observational data do not indicate any suppression in the 
matter power spectrum, roughly for the wavenumber k/h < 1 Mpc -1 . Hence, from 
the above expression, we can roughly obtain a constraint on the model parameter M 
as M > 10 eV. In the next subsection, we shall confirm this bound by numerically 
calculating the matter power spectrum. 

4.2 Numerical calculation 

4.2.1 Field description and background evolution 

In the numerical calculation, it will be practically difficult to use Eqs. ( p.6|) -(pT9|) 



because the sound speed of GDM increases beyond the speed of light as we go 
back in time. (Of course, this is not a physical problem but just a breakdown of 
the description.) Therefore, in the following, we shall solve X assuming a working 
Lagrangian to avoid this difficulty and follow the time evolution in a numerically 
stable way, instead of solving Eqs. (|3.6|)-(P^D. Let us consider a model in which 
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P(X) is given by 



P(X) = ^(X-Mr. (4.3) 



In this case M = M. The time evolution of X is given by 

dX _QX(X - M 4 ) 
d\na ~ 3X-M 4 ' 

and the energy density, pressure, and sound speed of GDM are given by 



(4.4) 



Pgdm = ^(3* 2 - 2XM 4 - M 8 ) , (4.5) 

^ g dm = ^(X-M 4 ) 2 , (4.6) 
2 _ X — M 4 

° s ~ 3X-M 4 ' (4 ' 7) 

respectively. The initial condition for the value of X(= X ) is fixed at present de- 
manding that GDM component should be accounting for the dark matter component 
observed today: 

°*- = rnm^^ - 2X ° M> - M8) = °- 3 (4 - 8) 

Then the evolution of X can be solved backward in time using Eq.( |4.4| ). 

As we have discussed in the previous section, we have to restrict our consideration 
to the regime \X — M 4 | /M 4 <C 1 and in this regime the sound speed evolves as 
c s oc yJ\X — M 4 | oc a -3 / 2 . Hence, if we look back the background evolution in 
AGDM universe, the condition \X — M 4 \/M 4 1 might be violated at some point 
in time. Let us define a critical scale factor when the condition \X — M 4 | /M 4 1 
is violated as 

|X(a cr ) - M 4 | /M 4 = 2c 2 (a cr ) = 0.1 . (4.9) 

From this equation, the critical scale factor a cr can be obtained in terms of M as 

4 ftt„ dm h 2 \ 1/3 A.0eV\ 4/3 . , 

aja K 1.0 x 10- 4 j . (4.10) 

Thus, it is found that for M > 1 eV the critical scale factor a cr is smaller than the 
scale factor at the matter-radiation equality, a eq . In the previous section, we have 
obtained the lower bound for M from the background evolution as M > 1 eV. Hence, 
if this bound is satisfied, we can use the low energy effective theory, at least, after 
the matter-radiation equality. Hereinafter, we consider the case for M > 1 eV. 



4 In Ref. p2| , the authors have also investigated cosmological implications for the dark matter 
model whose Lagrangian is given by Eq. (4.3), called "kinetic UDM" model in that paper. 
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This GDM description by X is equivalent to fluid-like description by Eqs. Q3.6|) - 
when X is near the condensate position, i.e., \X — M 4 \ <C M A . However, the 



sound speed will approach to c 2 = \ in this description and will never exceed the 



3 

r4i /n #4 



speed of light even if the energy density evolves to infinity as \X — M \ /M 3> 1, as 
is shown in Fig. [1|. The equation of state parameter, w = P g dm/ Pgdm, also approaches 
to 1/3. Therefore we can solve the system in a numerically stable way. It should be 
noted, however, that this description should be considered only as a regulator of the 
fluid description considered in the previous section, and it will be invalid when X is 
far away from the condensate position. Still one needs to interpret the results where 
a < a cr with caution. 

In fact, we have to extend our analysis before a cr because the density pertur- 
bations should be solved numerically from the early epoch when the corresponding 
Fourier modes are outside the cosmic horizon. Therefore, in the analysis in this sec- 
tion we extrapolate the behavior of GDM beyond a cr with a reasonable asuumption 
that the sound speed of GDM had been satulated as c 2 = 1/3 for a < a cr . This 
treatment corresponds to an assumption that the density perturbation of GDM can 
not grow in the radiation dominated era even when the low energy effective theory 
is broken. 

The corresponding Jeans wavenumber is plotted in Fig. ||. From the figure, we 
can see a turnover of the Jeans wavenumber kj for each M . This turnover of the 
Jeans wavenumber occurs at the transition of the sound speed of the GDM between 
c s = l/v3 and c s oc a~ 3 / 2 . The Jeans wavenumber is given by 

3 

^J,turn = ^^ a turn-f^( a turn) j (4-H) 

where we denote the scale factor at the turnover epoch by a turn . We can easily find 
that a turn ~ a cr , where a cr is given by Eq. (|4. 10Q . For M > 1 eV, of course, a tur n — «cr 
is smaller than a eq . 

4.2.2 Perturbation evolution 

Let us consider the evolution of the matter perturbations. In synchronous gauge, the 
growth of density contrast of GDM, 5 g d m , is described by the differential equations 
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<5gdm + (! + 3c 2 - 6w gdm ) "H<5g dm 

+ |^ 2 ( C s _ W gdm)(l - 6w gdm - 3w) 

+c 2 s k 2 + 3U{c 2 J - w' gdm ) j<5 gdm 



;i + w sdm )— - (1 + w gdm )(l - 3w gdm )n- - (4.12) 
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Figure 1: The evolution of the sound speed of GDM for various values of M. The 
blue short-dashed line is for cosmological constant plus ghost dark matter (AGDM) with 
M = 1CT 1 eV, the magenta dotted line is for the AGDM with M = 1.0 eV, the red 
solid line is for the AGDM with M = 10 eV and the cyan long-dashed line is for the 
AGDM with M = 10 2 eV. In all plots we set tt A = 0.7 and fi gdm = 0.3 (or Q cdm = 0.3 
for the ACDM). The vertical black dashed line represents the matter-radiation equality. 
The downward arrow shows the regime of the validity of the low energy effective theory 
(\X - M 4 | /M 4 < 0.1). 

h" + Uti = -3U 2 (l + c s 2 )5, (4.13) 

where a prime denotes a derivative with respect to conformal time, % = a' /a, h is 
the metric perturbation, 8, c s , w are total component's density fluctuation, sound 
speed and equation of state, respectively, denned as 

j = EM ffl= EA 3 = EA£i. ( 4.i4) 
Lp* l^Pi L pi 

The evolution of perturbations in the sound horizon depends on the time dependence 
of the sound speed. In the present model the sound speed is c 2 = 1/3 in the very 
early universe, and c 2 oc a -3 / 2 after some critical epoch. Therefore, the perturbation 
evolution in the sound horizon is divided into two cases. They are; case (a): w g d m = 
c 2 = 1/3, and case (b): w g( ± m ,c s 1. 

Note that the case (a) may not be a correct description of GDM scenario in 
the sense that the low energy effective theory is broken at large redshifts where 
a < a cr(^ Oeq). However, it would be natural to expect that the density perturbation 
in GDM can not grow very much if the sound speed of GDM had been saturated 
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Figure 2: The comoving Jeans wavenumber for various values of M as indicated in the fig- 
ure. From up to bottom, the lines are for M = 10, IkeV, 100, 10, 1 eV. The vertical dashed 
line shows the matter-radiation equality time. The turnover of the Jeans wavenumber oc- 
curs due to the change of the sound speed of GDM, which roughly shows the boundary of 
the regime of validity of the low energy effective theory. 

as c s ~ 0(1). This situation would be effectively taken into account by setting 
c 2 , = ~ because in this case the density perturbation can not grow at subhorizon 
scales. Therefore we shall assume case (a) in the following analysis to solve the 
perturbations in GDM for the energy scales in which the effective theory can not be 
applied. 

• case (a): w gdm = c 2 s = 1/3 

In this case, c' s ,w' dra = 0. So, differential equations(Eq. (|4.12| ), Eq. Q4.13|) ) are 



reduced to 

+ c 2 s k 2 5 gdm = --h", (4.15) 

h" + Hti = -m 2 5. (4.16) 

Roughly speaking, h" oc 1-L 2 5 from the second equation, and using the fact that 
% <C c s k in sound horizon, the first equation becomes 

^ dm + c^ 2 5 gdm = 0. (4.17) 

Therefore, the density perturbation of GDM in this case cannot grow and only os- 
cillate. 
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• case (b) : w g d m , c s C 1 

In this case, c' s = —^Hc 8 and w' gdm = SHwgdm- So, Eq. ( |4.12| ) becomes 

+ H5 Um + {fsk 2 ~ \u 2 ^ 5 gdm = 0. (4.18) 

Since H'(r}) <C c? s (r))k, we use WKB approximation 

<5 g dm = A(rf) exp J c s kdi^j , (4.19) 

to solve this equation. In this approximation, A'(rj) is roughly in proportion to % 
and c s k ^> H. So the term proportional to % 2 can be neglected. Then, the equation 
becomes 

{2c s A' + c s A + c s UA) ik exp ^ c s kdr^j = 0. (4.20) 

Then we obtain A oc a 1//4 since c' s = —^Hc s . Thus the growth rate of GDM for case 
(b) is 

5 oc a 1/4 , (4.21) 

in the sound horizon. Therefore, the density perturbation of GDM for case (b) 
gradually grows in the sound horizon. The above consideration can be generalized 
in the case where c Sj oc a~ n (n ^ 0). In this case, c' s = —n'Hc s , so the evolution of 
perturbation in sound horizon is 

5 oc a (n " 1)/2 , (4.22) 

which is derived from Eq. (|4.20|) . 

The evolution of GDM density perturbations in AGDM cosmology can be divided 
into four cases by comparing the wavelength of the modes and the Jeans length. Let 
us define aj(fc) as a scale factor when a mode becomes larger than the Jeans scale, 
namely, when the mode exits the sound horizon. We define the four cases from 
the smallest scale to the largest one, depending on aj(k) and the matter-radiation 
equality time denoted by follows; 

• case I (k < k JMrn ): 

the modes of interest never cross the sound horizon during the cosmic history, 

• case II (k JtCq > k > k Jtturn ): 

the modes of interest exit the sound horizon before the matter-radiation equality, 

• case III (do > a,j(k); k > kj, eq ): 

the modes of interest exit the sound horizon after the matter-radiation equality, 
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• case IV (aj(k) > a ; k > kj fiq ): 

the modes do not exit the sound horizon until today. 




Figure 3: Four cases (case I - IV) divided for the Jeans wavenumber for the AGDM with 
M = lOeV due to the difference of the evolution of the perturbation. Thick solid line shows 
the evolution of the Jeans wavenumber, vertical thin dashed line shows the matter-radiation 
equality. 

For example, Fig.|3] shows typical wavenumbers of four cases for the AGDM with 
M — 10 eV. If the mode belongs to the case I, the evolution of the mode is similar 
to that in the ordinary CDM one because the mode is always outside the sound 
horizon. However, if the modes belong to the cases II - IV, the modes enter the 
sound horizon at some epochs, where the perturbation evolution of GDM becomes 
different from that of CDM. Then, let us consider how the modes entering the sound 
horizon of GDM can be suppressed and estimate the fc-dependence of the matter 
power spectrum. In the M > leV models, the small scale perturbations cross the 
sound horizon in the radiation dominated era as shown in Fig|| However, the time 
of exiting the sound horizon dependends on the wavenumber of the mode k and also 
the cut-off scale M as mentioned above. 

From these considerations, one can estimate the power spectrum of GDM model. 
In the following, for simplicity, we neglect the phase of oscillations in the sound 
horizon. In fact, it is found to be important quantitatively and we will include 
this effect in the numerical calculation which will be presented later in order to 
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calculate evolution of perturbations and the matter power spectrum more precisely. 
We approximate <5 g d m oc a and oc log a outside the sound horizon in the matter and 
radiation dominated era (a cr < a < a eq ), respectively. Because of the fact that 
<5 g dm oc a 1 / 4 inside the sound horizon when a > a cr , we find 



' 4dm (oo) for fc < /cj iturn 

(aj/a cr ) 1/4 log (a C q/aj) (ao/a cq ) 5 gdm (a hc ) for fc Jjeq > k > k JMin 

(aj/a cr ) 1/4 (a /aj) <5 gdm (a hc ) for a > a 3 (k); k>k Jm 

^ (a /a cr ) 1/4 5 gdm (a hc ) for a 3 (k) > a ; fc > fcj, eq 



where we have introduced a nc as the scale factor when the modes of interest enter in 
the Hubble horizon. In these equations, only aj has k dependence, namely, aj oc k 2 
in the radiation era and aj oc k in the matter era. Meanwhile, in the standard cold 
dark matter cosmology, we find 5 c dm(«o) ~ log(a eq /ah c ) (ao/a e q) ^gdm(ahc)- Therefore, 
the k dependence of the GDM power spectrum at present can be found to be 



1 for k < k JMrn 



Pgdm.\k) 
Pcdm(k) 



(log k) for /c Jjeq > k> kj } turn / n a\ 

k~ 2 for a > aj(k); k > k Jm 
(log A;)" 2 for a^k) > a ; k > k Jm 



where we have used x~ 1 ^ i log(x) ~ 0(1) for x > 1. 

The analysis so far is based on the rough estimate and only appropriate for 
quantitative understanding. In order to evaluate the perturbation amplitude quanti- 
tatively and take into account the effects neglected in the above analysis, we calculate 
the evolution numerically using the modified CAMB code [23 . 

In Fig.f|, we depict the evolutions of density perturbations of GDM with different 
wavenumbers (k = 0.1, 1, 10, 100 Mpc" 1 ) for two different model parameters (M = 
20 eV and 100 eV). We also depict those of the baryon density for comparison. At 
largest scales (top panels) the evolution of GDM density perturbations is almost 
identical to that of standard CDM, because the mode is always outside the sound 
horizon. On the other hand, the evolutions of the modes at small scales are depicted 
in the third and fourth panels. In these panels the amplitude of GDM exhibits 
oscillations when the perturbations are inside the sound horizon. The growth of 
density perturbation of GDM inside the sound horizon for w g d m , c 2 s 1 is clearly 
seen, which has a dependence of a 1 / 4 as derived in our analytic estimate. We need, 
however, to see the second-top left panel (k — 1 Mpc -1 and M = 20 eV) with 
care. The panel shows the marginal case where the mode exits the sound horizon 
immediately after entering the horizon. In this case, we observe that the amplitude of 
5 g dm does not experience any logarithmic growth which is expected for the CDM case. 
We found that this brings a large difference in amplitudes at the present universe 
about an order of magnitude, which is seen by comparing with the second-top right 
panel. 
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Figure 4: The evolution of density perturbations of GDM and baryon fluid in Ghost 
condensation model. Left panels are for M = 20 eV model, right panels are for M = 100 eV 
model. The panels from top to bottom are for k = 0.1, 1, 10, 100 Mpc -1 . 



The net effect in the matter power spectrum is the deficit of power for k > kj tnTn , 
which is shown in Fig.|5]. From this figure, we can find that the matter power spectrum 
alters not at kj m as discussed in Sec. 4.1, but at fcj,t U m- That is because we have 
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extended the analysis beyond the critical epoch a cr in our numerical calculations, by 
extrapolating the behavior of GDM beyond a cr with a reasonable assumption that 
the sound speed of GDM had been satulated as c 2 s = 1/3. In this case the density 
perturbation of GDM only oscillates inside the sound horizon and does not experience 
any logarithmic growth during the radiation dominated era which is expected in the 
ACDM model. The power deficit at k > kj eq is more reasonable because it is a 
result derived mainly from the perturbation evolution at the subhorizon scales at the 
regime where the effective low-energy theory is applicable, i.e., (X — M 4 )/M 4 < 1 or 
a > a cr , under a simple assumption that the density perturbation of GDM can not 
grow deep in the radiation dominated era, even when the low energy effective theory 
is broken. 

The evolution of the baryon density perturbation is affected only through gravity. 
In the standard cosmology, it is known that after the fluctuation is erased by diffusion 
damping the baryon fluctuation should start to grow in time again with terminal 
velocity falling into the gravitational potential of clustered dark matter [p5[| . In 
the GDM model, however, the baryon fluid can not have terminal velocity because 
the dark matter density fluctuation oscillates in the sound horizon and hence the 
gravitational potential. As the GDM ceases oscillating, the amplitude of the density 
fluctuation of baryon fluid starts growing with oscillation, which is clearly seen in 
the left bottom panel. In any case, the baryon density fluctuation quickly catches 
up the dark matter density one after decoupling. 

The matter power spectrum of our universe can be probed through various cos- 
mological observations, such as clustering of galaxies ||26fl , cosmic shear [||, |7], ||], 



Lyman-a forest |3(J, [y], |32fl , and so on. At present, almost all observational data 
support the CDM paradigm, roughly for the wavenumber k/h < 1 Mpc -1 . In the 
models considered in this paper, this observational fact gives us a constraint on the 
model parameter M. The constraint is derived through the suppression of the mat- 
ter power spectrum, in the same way to obtain the constraints on the hot and/or 
warm dark matter models (or in other words, masses of neutrinos and/or warm dark 
matter particles) []53| . By looking at the matter power spectrum obtained in the 
present analysis (Fig.^j), we conclude that the model parameter should be 

M > 10 eV , (4.25) 



which is a stronger constraint than that obtained only from the background evolution 
considered in Sec. This result is consistent with the rough analytic estimate given 
by Eq. §3). 

Finally we should note that the way of suppression in the GDM models is different 
from that in hot or warm dark matter models. In the GDM models the suppression is 
the power law while for hot and warm models the power is exponentially suppressed. 
This fact might be used to distinguish between GDM and other dark matter models. 
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Figure 5: The matter power spectrum for GDM model. Red solid line is a matter power 
spectrum of ACDM model. The GDM power spectra are plotted for three models : M = 
1 eV (green dashed line), M = 20 eV (blue short-dashed line), M = 100 eV (pink dotted 
line). For modes k > A;j jturn (shown by a down arrow for M = 20 eV), the power spectrum 
of GDM model is suppressed compared to that of ACDM model. 




5. Summary and discussion 

We have investigated the possibility that ghost condensation may serve as an alterna- 
tive to dark matter, which had not been investigated in detail. In the present paper 
we have considered the AGDM universe, i.e. a late-time universe dominated by a 
cosmological constant and ghost dark matter. We have investigated the Friedmann- 
Robertson- Walker (FRW) background evolution and the large-scale structure of the 
AGDM universe, and have found a lower bound on the scale of spontaneous Lorentz 
breaking as M > 10 eV. This bound is compatible with the phenomenological upper 
bound M < 100 GeV known in the literature. 

As we have reviewed in Sec. |2|, ghost condensation is the simplest Higgs mech- 
anism for gravity in the sense that the number of Nambu-Goldstone boson is only 
one. The structure of the low-energy effective field theory is completely determined 
by the symmetry breaking pattern and this makes it possible for us to give robust 
predictions of the theory as far as the system is in the regime of validity of the effec- 
tive theory. In Sec. |2| we have also reviewed the infrared modification of linearized 
gravity as well as some non-linear dynamics and the phenomenological lower bound 
mentioned above. 

In Sec. |]we have provided a simple description of ghost dark matter and investi- 



gated the FRW background evolution. In subsection 3.1 we have shown that, under 



a certain condition, the background evolution and the behavior of large-scale pertur- 
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bations of ghost dark matter can be described by a fluid with the equation of state 
-Pgdm oc Pgdm/M 4 - This description has been used throughout this paper, except for 
numerical studies where the field picture turns out to be more convenient for some 
technical reasons. In subsection ^2, as a possible cosmological scenario relating the 
late time evolution of ghost dark matter to the early universe, we have considered 
a generation mechanism of ghost dark matter at the end of ghost inflation pi ]. In 



subsection |b\3| we have shown that the background FRW evolution in the AGDM 
universe is indistinguishable from that in the standard ACDM universe if M > 1 eV. 

In Sec.|] we have investigated the large-scale structure of the AGDM universe. 
Since the GDM has the effective sound speed unlike the standard CDM, small scale 
perturbations are suppressed. The suppression of the matter power spectrum occurs 
in a way different from that in models with hot dark matter or warm dark matter 
particles. We have given an analytic treatment to predict the matter power spec- 
trum observed today and have also calculated the power spectrum numerically. The 
analytic treatment can be extended for general dark matter models in which dark 
matter has a finite sound speed whose time dependence is the power law of scale 
factor. By comparing the GDM power spectrum with the observed one, one can get 
a lower bound on the scale M. The result we obtained is M > 10 eV. 

The constraint obtained in this paper can be improved further by observations 
of the matter power spectrum at smaller scales. For example, 21cm-line observations 
and/or the precise determination of the reionization epoch will provide us plenty of 
information about the matter power spectrum at smaller scales and hence stronger 
limits on the ghost condensation scale M. We leave these interesting subjects for 
future investigation. 

As suggested in JllJ , ghost condensation may provide an alternative explanation 
for the acceleration of the present universe if M ~ 10~ 3 eV. If the cosmological 
constant in the symmetric phase (with <p = 0) is zero then the effective cosmolog- 
ical constant in the broken phase, i.e. the ghost condensate, (with (f) = M 2 ) is 
positive 5 and of order 0(M 4 / Mpj) unless fine-tuned. Unfortunately, the condition 
M ~ 10 -3 eV is not compatible with the lower bound on M found in the present pa- 
per under the assumption that ghost dark matter is responsible for all dark matter in 
the universe. Therefore, it is not easy for the ghost condensate to be a simultaneous 
alternative to dark energy and dark matter unless fine-tuned. 

On the other hand, as shown in subsection \j.2[ ghost inflation is compatible with 
ghost dark matter. Ghost dark matter is naturally produced at the end of ghost 
inflation. Moreover, the lower bound ( [3.161 ) from ghost inflation can be satisfied 
simultaneously with not only the phenomenological upper bound ( |2.29|) but also the 
bound from ghost dark matter. Detailed investigation of the combination of ghost 
inflation and ghost dark matter is certainly worthwhile as a future work. 

5 Note that in the standard Higgs mechanism, the effective cosmological constant in the broken 
phase would be negative if the cosmological constant in the symmetric phase is zero. 
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